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TECHNICAL  PAPER 


SRM  PROPELLANT  AND  POLYMER  MATERIALS  STRUCTURAL  MODELING 

INTRODUCTION 


The  need  for  better  structural  characterization  and  analysis  of  the  solid  propellant  grain  is  dictated 
by  the  need  to  characterize  the  statics  and  dynamics  of  the  solid  rocket  boosters  as  a  system.  There,  pre¬ 
sently,  is  not  a  data  base  of  analysis  methods  and  supporting  test  data  which  allows  an  accurate  prediction 
of  the  propellant  or  bondline  strain  or  stress.  This  lack  of  test  and  analytical  structural  characterization  of 
propellants  results  in  unknown  safety  factors  and  a  general  inability  to  develop  criteria  for  the  quality 
control  necessary  for  manned  space  systems.  In  fact,  the  propellant  cumulative  structural  margins  for 
storage,  thermal  cycling,  transportation  induced  stress  cycling,  and  ignition  transients  have  never  been 
experimentally  verified  on  any  manned  solid  rocket  motor  (SRM). 

The  SRM  propellant,  insulation,  inhibitor,  liners,  and  seals  have  been  generally  characterized  as 
being  made  of  viscoelastic  materials.  Although  the  viscoelastic  classification  has  been  generally 
accepted,  close  examination  of  these  materials  reveal  that  they  are  either  more  complex  and  nonlinear 
than  classic  viscoelastic  models  or  the  actual  mechanisms  should  be  redefined  in  a  different  mathematical 
form.  The  following  investigation  reviews  and  evaluates  the  Space  Shuttle  SRM  propellant  structural 
analyses  utilizing  stress  relaxation  modulus  data.  The  stress  relaxation  test  data  is  examined  and  a  new 
math  model  is  proposed. 


STRESS  RELAXATION  MODULUS 


The  Stress  Relaxation  Modulus  is  obtained  by  inducing  a  constant  strain  or  shear  in  the  propellant 
test  specimen  and  measuring  the  reaction  boundary  forces  or  moments  with  respect  to  time.  The  stress 
relaxation  modulus  at  70°F  is  tabulated  versus  time  in  the  two  left  hand  columns  of  Table  1 .  The  data  of 
Table  1  was  obtained  from  Reference  4.  The  stress  relaxation  modulus  data  points,  for  time  greater  than 
10-®  min,  are  plotted  in  Figure  1 .  Table  2  and  Figure  2  show  the  temperature  dependent  time  shift  factor 
which  is  used  by  Reference  4  to  extrapolate  the  test  results  of  Figure  1  to  different  temperatures.  Time  is 
simply  divided  by  this  factor  which  has  the  effect  of  expanding  or  compressing  the  material  time  scale 
with  respect  to  the  actual  time  scale  of  the  transient  events.  The  modification  of  the  Stress  Relaxation 
Modulus  by  the  time  shift  parameter  is  consistent  with  the  assumption  fo  the  material  being  linearly 
viscoelastic.  Although  the  temperature  dependent  time  shift  parameter  is  generally  accepted,  the  author’s 
opinion  is  that  the  nonlinear  boundaries  of  the  real  material  must  be  defined  to  realistically  judge  the 
blanket  usage  of  this  assumption.  There  is  no  doubt  a  theoretical  region  where  the  temperature  dependent 
time  shift  parameter  is  suitable,  but  as  the  polymer  is  applied  theoretically  and  extended  to  nonlinear 
cases  other  than  stress  relaxation,  this  assumption  must  be  evaluated  closely. 


TABLE  1.  PRONY  SERIES  TERMS  FOR  SHEAR  MODULUS 


t 

(min.) 

G(t) 

psi 

Gi 

(psi) 

Log  \j 

10'^ 

80 

Go  =  75 

G,  =  8.24 

Log  X 1  =  6.3 

10^ 

92 

G:  =  14.51 

4.3 

10- 

130 

G3  =  29.4 

2.3 

1 

220 

G4  =  153.4 

0.3 

10-- 

420 

G5  =  231.3 

-1.7 

I0-* 

1,100 

Go  =  1,013.0 

-3.7 

10'^’ 

2,930 

G7  =  2,367.0 

-5.7 

10,000 

Gs  =  11,023. 

-7.7 

,0-10 

48,500 

Go  =  55,508 

-9.7 

10'- 

220,000 

G,o  =  247,120 

-11.7 

TABLE  2.  Log  aj  VERSUS  TEMPERATURE 

Log  aj 

Temperature,  °F 

-4.4 

145 

-2.0 

105 

0 

70 

2.1 

33 

4.4 

0 

6.7 

-30 

SHIFT  FACTORS  FROM 
LARGE  SPECIMEN  DATA 


50  100 

TEMPEKATUKM  (’’F) 

Figure  2.  Log  aj  versus  temperatui 
stress  relaxatio 


5 


Figure  3.  SRB  propellant  real  modulus  versus  frequency  and  dynamic  strains. 


Frequency  (coAi-),  Hz 

Figure  4.  Comparison  of  static  strain  levels  for  batch  No.  TP-Hl  148-9. 
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The  stress  relaxation  data  has  been  used  by  Reference  4  for  propellant  structural  analysis  of  the 
Space  Shuttle  ARM.  One  method  is  simply  for  the  engineer  to  assume  or  select  a  time  during  an  event 
such  as  the  long  term  storage,  then  the  modulus  is  defined  by  the  stress  relaxation  data  and  a  linear  elastic 
analysis  is  completed  for  that  one  instant  in  time.  It  should  be  remembered  that  the  stress  relaxation  mod¬ 
ulus  is  only  valid  for  constant  strain  while  many  actual  problems  are  really  transient  in  nature.  It  has  been 
proven  in  test  that  under  transient  conditions  the  modulus  is  strain  rate  and  static  strain  dependent  (Figs.  3 
and  4). 


The  second  method  of  analysis  is  to  use  a  finite  element  linear  viscoelastic  transient  computer 
code.  The  code  currently  used  by  Reference  4  uses  a  Prony  series  to  model  the  modulus  of  the  assumed 
linear  viscoelastic  material,  as  shown  below  (Table  1): 


=  Go  +  G,  e-‘^^'  +  G2  +  G3  e-‘^^3  -t- 


G(t) 


Reference  4  used  only  the  stress  relaxation  data  to  define  the  terms  in  this  series,  and  thus,  again  the 
model  as  currently  used  may  only  be  valid  for  the  constant  strain  condition,  which  is  hardly  a  typical 
transient  condition. 

The  Prony  series  material  model  (Fig.  5)  can  be  defined  conceptually  by  a  very  simple  schematic 
consisting  of  a  spring  in  parallel  with  a  number  of  Maxwell  elements,  as  in  Reference  I .  The  Maxwell 
elements  are  simply  a  spring  in  series  with  a  viscous  dash  pod  damper.  This  is  the  simplest  and  most 
numerically  convenient  viscoelastic  model  for  a  material.  A  time  dependent  analysis  is  then  done  with 
this  model  assuming  linear  superposition  and  ignoring  a  variety  of  factors  which  can  cause  the  effective 
modulus  to  shift  by  factors  of  two  and  more. 


Figure  5.  Conceptual  schematic  of  Prony  Series  material  model. 
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INHERENT  PRONY  SERIES  SRM  MATH  MODEL  WEAKNESSES 


The  Space  Shuttle  SRM  Prony  Series  structural  model  has  several  inherent  weaknesses  numer¬ 
ically,  experimentally,  and  theoretically. 

1 .  A  large  number  of  terms  are  used  by  Reference  4  to  match  a  comparable  number  of  stress 
relaxation  test  points.  (11  terms  with  a  total  of  21  constants  as  shown  in  Table  1.) 

2.  The  series  is  extremely  ill  conditioned.  There  are  many  combinations  of  G| ,  G2,  G3,  . . . ,  and 
X.1,  X2>  ^3,  ■■■■>  that  can  be  fit  to  data  generated  by  exact  functions  to  four  significant  figures.  The  intro¬ 
duction  of  experimental  noise  makes  the  operation  hopeless  [2]. 

3.  The  series,  as  fit  in  Reference  4,  does  not  have  a  smooth  uniform  curve  between  the  stress 
relaxation  data  points  even  though  it  passes  through  the  data  points  (Fig.  6). 

4.  The  derivation  of  the  dynamic  modulus  from  the  Prony  Series  using  the  assumption  of  linear 
superposition  does  not  match  experimental  data  in  magnitude  or  trend  (Figs.  7  and  8).  The  real  part  of  the 
complex  modulus  is  labeled  as  G'.  The  imaginary  part  of  the  complex  modulus  is  labeled  as  G".  The 
magnitude  of  the  complex  modulus  is  labeled  as  G*. 

5.  Current  linear  viscoelastic  codes  and  Prony  Series  models  do  not  account  for  the  effects  of 
static  strain  and  dynamic  strain  (strain  rate)  on  dynamic  modulus  (Figs.  3  and  4). 

6.  The  original  Space  Shuttle  SRM  propellant  model  was  based  on  test  data  only  valid  for  con¬ 
stant  strain  (Reference  4). 

7.  The  series  fit  to  constant  strain  test  data  does  not  have  in  itself  any  relation  to  the  material 
response  unless  >12,  ^3,  •••  are  selected  properly  based  on  understanding  of  the  response  phenomena 
and  the  material  mechanisms. 

8.  The  region  of  the  lowest  propellant  safety  factor  is  the  forwrd  segment  star  pattern.  Although 
the  constant  strain  test  assumption  may  be  very  good  for  the  round  smooth  bore  segments  under  constant 
pressure,  the  assumption  is  not  valid  for  the  star  pattern  region. 
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Figure  6.  Shear  modulus  G(t)  versus  log  time,  Prony  Series  fit. 


SRM  REDESIGN  ANALYSIS  TECHNIQUES 


The  Space  Shuttle  SRM  redesign  was  checked  by  an  independent  and  direct  use  of  propellant  test 
data.  The  redesign  structural  analysts  have  modeled  the  propellants  and  insulation  as  elastic  isotropic 
materials  with  an  analysis  provision  or  capability  of  large  strains  and  a  variational  stress  relation  which 
allows  more  accurate  analyses  of  nearly  incompressible  materials  (References  10  and  12).  The  nonlinear 
variation  in  the  material  properties  was  accounted  for  by  reviewing  test  data  and  selecting  isotropic 
properties  which  were  the  best  approximation  for  the  particular  load  case.  Proper  review  of  the  validity  of 
these  analyses  requires  a  presentation  of  the  utilized  test  data  and  a  review  of  the  selection  of  the  equiv¬ 
alent  elastic  isotropic  material  data  (Reference  10).  There  was  no  presentation  of  a  test  verification  of  the 
material  property  selection  and  the  resulting  finite  element  math  model. 


The  analysis  of  the  operational  load  cases  involving  a  combination  of  load  cases  of  a  different 
nature,  such  as  pressure,  thermal,  and  acceleration,  are  analyzed  separately  with  totally  different  material 
properties  selected  for  each  case.  These  separate  cases  are  then  linearly  added  for  the  operational  condi¬ 
tion.  The  linear  superposition  with  the  limiting  constraints  and  assumptions  of  the  above  analysis  is  con¬ 
sistent  for  small  strains.  It  is  not  mathematically  correct  to  apply  linear  superposition  to  a  large  strain 
analysis.  There  has  been  no  test  verification  of  the  actual  superposition  relations  of  these  nonlinear 
materials. 

The  dynamic  transportation  analysis  was  termed  very  crude  by  the  author  of  the  reviewed  report 
(Reference  10)  and  in  his  own  words  should  be  considered  only  an  approximate  with  a  more  detailed 
analysis  recommended.  This  analysis  first  calculated  a  frequency  of  the  region  in  question  by  iterating 
with  assumed  shear  modulus  and  dynamic  modulus  test  data  from  Reference  8.  It  should  be  cautioned 
that  due  to  the  sensitivity  of  the  dynamic  modulus  to  static  strain  and  strain  rate  this  test  data  is  probably 
only  valid  for  the  test  conditions  and  as  documented  can  not  be  accurately  extended  within  a  factor  of  4  to 
the  SRM.  After  a  frequency  was  found,  a  g  level  or  environmental  specification  was  obtained  from 
Reference  1 1 .  A  load  and  resulting  strains  were  calculated  using  a  linear  elastic  static  solution  with 
factors  applied  for  transmissibility  at  resonance  with  a  33  percent  damping.  Although  the  dynamic  mod¬ 
ulus  test  data  can  be  questioned,  a  state  of  the  art  solution  using  this  test  data  would  have  been  a  frequ¬ 
ency  response  analysis  using  the  transportation  environment  as  driving  forces. 

The  above  analysis  reflects  the  need  for  the  development  of  a  test  verifiable  math  model  for  the 
SRM  propellant  and  related  polymer  materials.  The  above  reviewed  work  needs  a  demonstrated  and 
documented  test  verification.  Also,  the  analyses  is  totally  dependent  on  the  engineer’s  selection  of 
material  data  with  no  test  verified  and  NASA  endorsed  standards. 
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NEW  STRESS  RELAXATION  MODEL 


The  weaknesses  of  the  fit  experimental  data  to  the  Prony  Series  prompted  an  investigation  of  the 
test  data  as  shown  in  Table  1 .  The  shear  modulus  versus  time  of  Figure  1  is  replotted  in  Figure  9  on  the 
log  log  scale.  It  is  seen  that  the  table  data  points  can  be  fit  with  a  single  parabolic  curve.  Both  Figures  1 
and  9  show  that  the  shear  modulus  is  approaching  a  constant  static  value.  A  further  simplification  can  be 
made  by  subtracting  out  the  static  modulus  and  replotting  the  dynamic  shear  modulus  versus  time  on  a 
log  log  scale,  as  shown  in  Figure  10.  As  plotted  in  Figure  10,  the  stress  relaxation  data  looks  like  it  can  be 
fit  with  straight  lines.  A  straight  line  would  do  a  good  job  of  fitting  data  taken  at  times  greater  than  10'^ 
min.  The  data  points  for  less  than  10'^  min  either  indicate  a  different  phenomena  or  more  likely  problems 
with  testing  or  data  reduction  techniques  employed  in  Figure  1  and  Table  1  from  Reference  4.  The  stress 
relaxation  modulus  can  be  related  by  the  following  simple  equation. 


G(t)  =  Go  +  G,  r" 


where 


t  =  time  in  seconds 
Go  =  75  psi 
G|  =  331  psi 
n  =  0.2255  . 


Figure  1 1  is  a  plot  of  the  above  equation  and  should  be  compared  to  Figure  6.  The  above  equation  can 
stilt  be  used  with  the  temperature  dependent  time  shift  parameter. 

The  theoretical  dynamic  shear  modulus  can  be  derived  from  the  above  power  equation,  assuming 
linear  superposition  as  follows. 


Real  Modulus 


-JC, 

/ 


—  0)  J  sin(o)s)  G(s)  ds  [Ref.  1] 
o 


X 

!,(»)/[ 


G'(a))  =  Go  +  G|  fa)  J  [sin  (fa)s)/s"]  ds 

o 


G'(fa))  =  Go  +  G|  fa)"  r(l— n)  cos(mr/2) 
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Figure  9.  Log  shear  modulus  G(t)  versus  log  time. 


1000000 


Imaginary  Modulus 


G"((o)  = 


cos  (cos)  G(s)  ds 


[Ref.  1] 


/ 

o 


G"(a))  =  G,  to  I  [cos  ((os)/s")  ds 


G"((o)  =  G|  to"  r(l-n)  sin(mr/2) 


Complex  Modulus 


G*((o)  =  "v/C'Cco)^  +  G"(to)^ 


The  above  result  is  plotted  on  Figure  12. 

The  above  derived  dynamic  modulus  can  be  compared  to  the  theoretically  defined  dynamic 
modulus  based  on  the  Prony  Series  [1]. 


Gj  co^  X-i^ 
1  +  (0^ 


Prony  Real  Modulus 

G'(w)  =  Go  + 


N 

S 

i=  1 


Prony  Imaginary  Modulus 


N 

Gj  CO  X-i 


Prony  Complex  Modulus 

G*((o)  =  Vg'((o)2  +  G"(a))^ 

The  above  result  was  plotted  on  Figure  8. 
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Although  the  power  equation  theoretically  derived  dynamic  modulus  (Fig.  12)  was  based  on  the 
same  stress  relaxation  data,  the  results  are  very  different  from  the  result  shown  m  Figure  8.  The  Figure 
results  were  derived  from  the  Prony  Series.  It  can  be  seen  that  the  derived  dynamic  shear  modulus  plotted 
in  Figures  8  and  12  do  not  compare  to  the  experimental  results  of  Figure  7  very  well.  The  experimental 
test  data  obtained  from  Reference  8  as  shown  in  Figure  7  was  documented  as  being  driven  with  an 
electromechanical  vibrator  at  a  constant  displacement.  If  the  experimental  results  of  Figure  7  were 
actually  the  result  of  being  driven  at  a  constant  maximum  displacement  then  as  the  frequency  increased 
the  velocity  and  strain  rate  increased  in  proportion.  Figures  3  and  4  show  that  the  dynamic  shear  modulus 
is  significantly  shifted  by  changes  in  the  dynamic  and  static  strains.  If  the  experimental  data  had  been 
driven  at  a  constant  maximum  velocity  and  strain  rate,  the  experimental  results  probably  would  have 

been  different. 

The  theoretically  derived  dynamic  shear  modulus  based  on  either  the  Prony  Series  or  the  Power 
Equation,  as  they  currently  stand  cannot  predict  this  shift  with  dynamic  and  static  strain.  These  effects 
must  be  included  and  any  testing  or  analysis  done  without  the  proper  attention  to  these  parameters  is  o 
questionable  nature.  If  the  stress  relaxation  modulus  is  to  be  related  to  creep  tests,  dynamic  modulus 
tests,  and  strain  rate  tests,  these  two  parameters  must  be  included  and  checked  for  temperature  depen 
ency  and  validity  of  the  current  time  shift  parameter. 

The  author  would  like  to  propose  a  new  model  for  the  material  properties  of  the  propellant  and 
other  polymer  materials.  The  modulus  can  be  defined  by  a  summation  of  terms  where  each  te™  ^ 
function  of  multipliers.  Each  multiplier  is  of  two  possible  forms.  The  first  is  a  power  equation  and  the 
second  is  an  exponential  of  a  power  equation. 

An  example  of  the  series  can  be 


G  = 


Gstatic 


+  G 


Dynamic 


where,  Gstatic  reduce  to  a  constant  such  as  75  psi  for  the  propellant  and  Goynamic  inay  be  defined  as  a 
function  of  time,  temperature,  static  strain,  dynamic  strain,  and  other  terms  as  needed. 


GDy„amic(t,T,e,e)  =  M(t)  M(T)  M(8)  M(8) 


where  M(X)  takes  two  possible  forms: 


M(X)  =  C,  MX) 
or 

M(X)  =  C, 
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where 


Xq  may  be  equal  to  X  or  an  original  value  of  X 

Xr  is  a  X  reference  value 

Cn.  is  a  constant  in  the  numerator 

Cd,  is  a  constant  in  the  denominator 

Px  is  the  power  of  the  function  of  X  . 


As  an  example,  the  first  propellant  model  would  probably  use  a  p  function  of  order  one. 

tX)  =  +  Ck,  (X-Xr) 

^  Cd„  +  Cd,  (Xo  -  Xr) 


The  first  choice  for  the  propellant  multiplier  would  be; 

M(X)  =  C,  (B,(X) 

Applying  the  model  to  the  current  test  results  for  stress  relaxation  modulus 
M(t)  =  331  t-°-225 

|-  -]  8,5275 

I  285 

~  T  +  220  (adapted  from  Fig.  2) 
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M(e)  not  yet  defined 


Mfe)  not  yet  defined 


The  exponential  form  of  the  multiplier  could  model  the  Prony  Series  form  and  also  adapt  to  past  accepted 
models  of  different  polymers.  As  an  example,  the  temperature  dependent  time  shift  parameter  that  has 
been  accepted  for  polyisobutylene,  natural  rubber,  polyurethane  elastomer,  polystyrene,  and  poly(ethyl 
methacrylate)  is  as  follows  from  Reference  1 . 


log  Ctj  = 


-C3(T-Tr) 
C4  +  T-TR 


The  temperature  dependent  multiplier  would  be  of  the  following  form  for  this  relation. 


CONCLUSIONS 


A  new  math  model  for  stress  relaxation  has  been  proposed.  This  model  has  been  extended  and 
generalized  into  a  form  utilizing  multipliers  with  each  multiplier  incorporating  the  contribution  of  one 
important  parameter  such  as  time,  temperature,  strain  rate,  static  strain,  etc.  The  final  definition  of  all  the 
necessary  terms  and  the  general  application  of  this  model  must  be  developed  and  verified  by  detailed 
consistent  test  program  for  the  materials  in  question  as  in  Reference  9.  Although,  it  is  hoped  this  model 
will  be  very  successful,  it  is  merely  the  next  step  which  is  highly  dependent  on  practical  test  results,  and 
it  must  be  suitably  test  verified. 
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APPENDIX 


Additional  plots  used  in  the  evaluation  of  the  equation  forms  of  the  stress  relaxation  modulus  and 
the  dynamic  modulus  are  given  in  Figures  13  through  23. 
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